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Abstract 

As a prototype of an evolution equation we consider the Schrodinger equation 

i(d/dt)'B(t)=H'B(t), H = H 0 + V(x) 

for the Hilbert space valued function 'k(-) : M —> TL = L‘ 2 (M u ) which describes 
the state of the system at time t in space dimension v > 2. The kinetic energy 
operator Hq may be Hq = — (l/2m) A (nonrelativistic quantum mechanics for 
a particle of mass m), Ho = \J —A + m 2 (relativistic kinematics, Klein-Gordon 
equation), the Dirac operator, or ..., while the potential V (x) —> 0 suitably as 
|x| —> oo. 

We present a geometrical approach to the inverse scattering problem. For 
given scattering operator S we show uniqueness of the potential, we give ex¬ 
plicit limits of the high-energy behavior of the scattering operator, and we give 
reconstruction formulas for the potential. 

Our mathematical proofs closely follow physical intuition. A key observation 
is that at high energies translation of wave packets dominates over spreading 
during the interaction time. Extensions of the method cover e.g. Schrodinger 
operators with magnetic fields, multiparticle systems, and wave equations. 


1 Introduction, the Schrodinger Equation 

The Schrodinger equation is a linear evolution equation for a function of time 1 6 R 
with values in a state space (phase space) TL which is a Hilbert space: 

T(-) : R -> TL. 

The initial value problem reads 

= T(0) = T, (1.1) 

with a linear operator H acting on TL. (We use units where Planck’s constant h = 1.) 
This type of equation includes as special cases nonrelativistic and relativistic quan¬ 
tum mechanics, the Dirac equation, the linear wave equation (with the usual method 
to transform a second order equation into a first order system), and other evolution 

x To appear in: Understanding Physics, A.K. Richter ed., Copernicus Gesellschaft, Katlenburg- 
Lindau, 1998, pp. 31-; ISBN 3-9804862-2-2 (Proceedings Bonn 1996). 
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equations. Splitting off the factor i is just for convenience because in important appli¬ 
cations H is symmetric. Since the operator H is typically unbounded, care is needed 
to ensure solvability of the equation. In the models mentioned above the operator H is 
self-adjoint on a suitably chosen domain T>(H). Then Stone’s theorem (or the spectral 
theorem and functional calculus, see e.g. ]16| ) ensure that the exponential exp {—itH} 
is a well defined unitary operator for all i 6 K and that 

T(t) = e~ itH ^ (1.2) 

is the unique global solution of the initial value problem ra¬ 
in the following we will describe our geometrical approach to the inverse problem 
for the Schrodinger equation as an equation which describes the motion of particles 
according to the laws of quantum mechanics. We will exploit physical intuition to 
help us solve the mathematical problems. The methods and results carry over to other 
evolution equations as well. 

We will discuss the differences of nonrelativistic and relativistic kinematics as far 
as they are relevant here. The time scales for interaction and for spreading of wave 
functions differ at high energies. This implies the simplicity of the leading behavior 
of the scattering operator because only the translational part of the time evolution 
matters as long as the interaction is strong. We obtain explicit formulas for the high 
energy scattering operator which can be used to reconstruct the potential uniquely. In 
the present paper we want to explain why the statements are true and how physical 
intuition and mathematical proofs are closely analogous. While we give all major steps 
of the proofs for two typical examples we refer to the papers for some more technical 
estimates and further examples. 

Acknowledgement. This paper is dedicated to Wolfgang Kundt who has strongly 
influenced my view of science. Part of the material was prepared during my stay at 
the Institute for Advanced Study, Pinceton, NJ, USA. I gratefully acknowledge the 
hospitality and support. 

2 Particles in Quantum Mechanics 

We describe the state of a quantum mechanical particle in u -dimensional space by a 
normalized vector T G hi. (Due to the superposition principle we may restrict ourselves 
to pure states and choose a normalized vector instead of the equivalence class of vectors 
differing by a global phase factor.) The vector can conveniently be represented in various 
ways (similar to basis changes in linear algebra), e.g. by a square integrable function 
ijj(-) G L 2 (IRP, dx) with volume measure dx. Instead of a function depending on the 
position or configuration space variable xel" one can use its Fourier transform 

$(■) G L 2 (W, dp), V>(p) := (2t r)-^ 2 J dx e" ipx ^(x) (2.1) 

depending on the momentum variable p G BP with normalization 

||T|| 2 = j dx |^(x)| 2 = j dp |^(p)| 2 = 1. 

We use for the abstract state vector a capital letter T, for its representation as a 
function of position -0(x), or its momentum space wave function p), respectively, 
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and write 

H <—* L 2 (W, dx ) ^—♦ L 2 (R u , dp) 

T *—> ip (x) ♦—> ^(p) (2.2) 

to indicate the switching between representations. 

For a given state T the probability measures /j x on conhguration space and /i p 
on momentum space, respectively, 

fi x (A) = f dx |^(x)| 2 and /Jtp(B) = [ dp \ip(p )\ 2 (2.3) 

J A J B 

describe the probabilities to find the particle in the (Lebesgue measurable) subset A C 
R" of conhguration space or 5 C 1" of momentum space. One may visualize such 
a state as a cloud of very many particles where /i x (A) describes the fraction of them 
which have their position in A and, similarly, /i p (B) is the fraction with momentum 
in B. Such a state is also called a wave packet. Classical point particles with a 5-like 
distribution in position and/or momentum space are impossible in quantum mechanics. 
The state space does not contain such idealized objects, in agreement with observations. 

We extend the triple of representations of state vectors to the linear operators acting 
on them. 

Abstract opera- Action on conhguration Action on momentum 

tor acting on 7i space wave functions space wave functions 

The Fourier transformation ( |2.1|) interchanges differentiation and multiplication of a 
function with its argument. Thus we obtain for the position and momentum operators, 
respectively, 

x <—» x *—» ?'V P , (2.4) 

p <—* —zV x <—> p. (2.5) 

In our notation sometimes we do not distinguish between the abstract operator on 7i 
and its representaion as a multiplication operator on the corresponding space on which 
it is “diagonal”. 

If the forces acting on the particle are described as the negative gradient of a po¬ 
tential function V(x) (conservative mechanical system) then the generator H of the 
time evolution, the Hamiltonian or Schrodinger operator , is the energy operator 

H = H 0 + V(x) (2.6) 

which is a sum of the kinetic energy operator Ho - responsible for the kinematics - and 
the real valued potential energy which determines the dynamics. The potential should 
decrease suitably as |x| —> oo, see Section [| for precise conditions on V. 

3 Kinematics 

The kinetic energy operator or free Hamiltonian H 0 usually is a function hfo(p) °f the 
momentum of the particle. We will study two typical cases, nonrelativistic (NR) and 
relativistic (Rel) kinematics. In the Erst case 
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It acts as a multiplication operator on 0 and as a differential operator on 0: 

H 0 $ < -+ (H 0 &)(x) = ~^-(A0)(x) ♦-+ H o (p)0(p) = ^p 2 0(p). 

Generally, the velocity operator is the change of position in time: 

= i[H 0 , x] = V p tfo(p) , (3.2) 

a function of the momentum operator. In the nonrelativistic case it is 


v(p) =-e itHo xe- itHo 
dt 


NR: v (p) = — . (3.3) 

m 

Note that the velocity is unbounded, the speed tends to infinity if the kinetic energy or 
the momentum does so. 

Let us now turn to the relativistic case corresponding e.g. to the Klein-Gordon 
equation. 

Rel: H 0 ( p) = \J p 2 c 2 + m 2 c 4 = \/p 2 + m 2 (3.4) 


if we use units of measurement such that the speed of light c = 1. Here the velocity 
operator is 


Rel: v(p) = VpiLo(p) = c 


Pc 


P 


\J p 2 c 2 T m 2 c 4 ^/p 2 + m 2 


(3.5) 


In this case the speed is bounded by c = 1 (the speed of light). 

The free time evolution operator is a simple multiplication operator in momentum 
space [and a complicated oscillating convolution operator in configuration space] 

e - itHo $ <—> (e- itHo 0)(x) <—> e-^ o(p) 0(p). (3.6) 


While for short times the free classical and quantum time evolutions differ considerably 
they behave similarly for large times. Asymptotically, the distribution in configuration 
space of a quantum state is in good approximation the same as that of the corresponding 
cloud of free classical particles, of the “classical wave packet”. For later applications 
we study a particular family of states <hp with compact momentum support around a 
very large average momentum peT. The operator exp(ipx), a unitary function of 
the position operator x, shifts a state in momentum space by p: 

$0 <— 0o(') <— 0o(") G CHM 17 ) (3.7) 

4>p =e* px <h 0 ♦—+ 0p(x) = e* px 0 o (x) «—+ 0 p (p) = 0 o (p - p). (3.8) 

Since 0 O (-) G «S(M I/ ), the Schwartz space of rapidly decreasing functions, these states are 
well localized in configuration space, too, uniformly in p. They have average velocities 
around v(p) G M v , where 


p/m 


NR, 


v(p) = Vi7 0 (p) =: i>(p) w 


p/\/p 2 + m 2 Rel. 


(3.9) 
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In the given examples the unit direction vector lj = v(p)/|v(p)| = p/|p| of the velocity 
v(p) has the same direction as p and the speed u(p) = |v(p)| depends on |p| only. 

Consider in the following only “large” p such that the minimal velocity in the 
support of (ftp is at least 2/3 rd of the average velocity: 


inf 



p G supp 




(3.10) 


In our examples e.g. |p| > 3sup{|p| | p G supp0o} will do. As we are mainly 
interested in the high-energy behavior |p| —> oo this restriction ( |3.10| ) is harmless, it 
excludes components in the state with particles of low or zero velocities which would 
require special treatment below. 

In our context we have to control the localization in configuration space of freely 
evolving wave packets. This depends mainly on the support of the state in velocity 
(momentum) space. Therefore, we have chosen compactly supported momentum space 
wave functions. Then in configuration space the states cannot have compact support as 
well but rapid falloff is sufficient there. A special case of such propagation properties 
(or non-propagation) of quantum wave packets for long times | 


17 


is 


J dx ( e -' !W O0_)(x)| 2 < 

x|<t t)(p)/2 


const($ 0 , n) 
(i + |t»(p)|)” 


(3.11) 


for any n G N uniformly for large p. A classical free particle which starts at time 
0 from the origin and has momentum p G supp (ftp will be localized at time t in the 
region 


x(t) e I x = 1 v (p) 


p G suppi 


C {x 


fv(p)| < u(p)/3}. 


(3.12) 


The “classically forbidden” region |x| < tv( p)/2 is separated from the “allowed region” 
by at least tv( p)/6. The state mainly propagates within the classically allowed region 
which moves away from the origin with a positive minimal speed. The “quantum 
tails” of the wave packet in the classically forbidden region do not vanish, nevertheless, 
they decay very fast in time, both in the future and past. This is physically and 
mathematically in close analogy to rays versus waves in optics. While the shadow 
behind an obstacle is not totally black due to diffraction it is, nevertheless, quite dark 
away from the region which can be reached by straight rays (the role of the increasing 
separation tv( p)/6 ). We will need below only the estimate of non-propagation (|3 .1 1| ) . 
It is proved with a stationary phase estimate of an integral with rapidly oscillating 
integrand. 

Kinematics has a strong influence on the high-energy behavior of the scattering 
strength. Wave packets will be influenced by the potential only as long as they are 
mainly localized in the region where the potential is strong. Due to the propagation 
property this is essentially a finite time interval, the interaction time, proportional to 
u(p) -1 . Thus the strength of the interaction (which is proportional to the duration of 
the interaction) will vanish as |p| —» oo in the nonrelativistic case and it will have a 
finite limit in the relativistic case. This is a crucial difference of the two cases. 













V. Enss 


6 


Inverse Scattering 


4 Dynamics 


Let us now turn to the interacting (perturbed) time evolution 
e~ itH tf, H = H 0 + V(x). 


(4.1) 


The free Hamiltonian H 0 is self-adjoint on its domain 

V(H 0 ) = (T £ H I 0(p) and H 0 ( p) 0(p) £ L 2 (M y , dp)}. 


The configuration space wave functions lie in the Sobolev spaces hL 2,2 (M y ) in the 
nonrelativistic case and in H /1,2 (M 1 ') in the relativistic case, respectively. First we will 
consider short-range potentials V(x) = Y s (x) which are roughly those which decrease 
at least like |x|”( 1+£ ), e > 0, as |x| —> oo. More precisely, the set of short-range 
potentials is 


V s 



sup \V S 

|x| >R 


dR < 


oo 


(4.2) 


For simplicity of presentation we will restrict ourselves in this paper to bounded poten¬ 
tials. The sum is self-adjoint on the domain T>(H) = T>(Hq). Singular potentials can 
easily be included using standard techniques ||16|| . This covers the physically relevant 
local singularities like l/|x| for Coulomb or Yukawa potentials. One simply has to 
regularize the singularities by a free resolvent or by a function /(p) of the momentum 
operator which has compact support and is the identity on <3>o • In addition, with suit¬ 
able adjustments we will treat long-range potentials V e which decrease towards infinity 
e.g. as slowly as the Coulomb potential, see Section [U]. Then 

V (x) = Y s (x) + I/ £ (x), 1Y e V s . (4.3) 


In the present context a short-range potential behaves similarly to a compactly 
supported one. Depending on the required accuracy it is essentially concentrated in a 
ball of some radius R around the origin. 

The influence on the particle by the force — W (x) is relevant only as long as 
the particle is essentially localized in the interaction region, i.e. where the potential is 
strong. Correspondingly, the state space can be split into two orthogonal components. 
The bound states remain forever under the influence of the potential, they constitute 
the pure point spectral subspace which is spanned by the eigenvectors of H. These 
states remain localized uniformly in time. Orthogonal to these are the scattering states 
which form the continuous spectral subspace H cont (H) = (eigenvectors of H}- 1 -, they 
leave the interaction region for large times (see e.g. Theorem XI. 115 in |L7|). The latter 
are studied in scattering theory. 


5 Scattering 

As a general reference for mathematical scattering theory see e.g. 0 or m for the 
Dirac equation. It is well known that for short-range potentials (H = Ho + V s ) the 
asymptotic motion of scattering states is an essentially free motion: For any scattering 
state T £ 7i cont (H) there exist free asymptotic configurations < L ± £ 7i such that 

\\ e -it[H Q +V s ] _ e ~itH 0 ([ji 


0 as t —> ±oo. 


(5.1) 
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This is usually called asymptotic completeness of the wave operators. Similarly, for 
any incoming configuration or outgoing <3> + there is a corresponding state T G 
Ti. cont (H) such that (|5.1|) holds (existence of wave operators). These statements can 
be proved for any dimension v of configuration space using propagation estimates like 
(|3.11|) . The proof of ( [7.2|) below is similar. 

A convenient tool to describe scattering is the scattering operator S which maps 
an incoming conhguration 4> _ to the corresponding outgoing conhguration <h + of the 
same state T. For given $' let 


T= lim e * f - [Ho+ys] e-"- 1 * 0 4>" and 

t_ —> —OO 

$+ = lim e it + H ° e - it +l H ° +ys ] T. 

t |_ —>OG 


Then 


£ ) • _ e it+H 0 e ~it + [H 0 +V s ] e it_[H 0 +V s ] £ ~it_ Hq 

S s-lim satisfies S <3> _ = <h + . 

t |_ —>oo 
t_ —>—oo 

For microscopic particles for which quantum mechanics is an adequate description one 
cannot really observe more details of the scattering process than those encoded in the 
scattering operator. We denote the mapping 

V S ^L(H), V S ^S=S(V S ) (5.4) 

as the scattering map from short-range potentials to bounded (unitary) scattering op¬ 
erators on the Hilbert space of asymptotic configurations. 

The direct problem of scattering theory is to determine for a given potential V 
the scattering operator while the inverse problem is to determine the potential(s) if 
the scattering operator or part of it is known. 


(5.2) 

(5.3) 


6 Uniqueness of the Potential 


We denote by F(H 0 > E ) the multiplication operator in momentum space with the 
characteristic function of the set {p G | H 0 ( p) > E}, i.e. the spectral projection of 
the kinetic energy operator to energies above E. The main results about uniqueness are 
of the following form. They are a corollary of the asymptotic behavior of the scattering 
operator shown below. 

Theorem 6.1 The scattering map S: V s —> L(Tt) is injective. Actually, the high- 
energy part of the scattering operator alone: S F(H 0 > E), E arbitrarily large, deter¬ 
mines the short-range potential uniquely. 


In the nonrelativistic case similar results go back to Faddeev (1956) and Berezanski 
(1958), the strongest results were by Salto (1984), all using time-independent methods. 
Our geometrical time-dependent proof covers a wider class of potentials and, more 
importantly, it is simpler. It has been extended to show analogous results for long- 
range potentials, magnetic fields, multiparticle systems, for the Dirac equation, and 
other systems, see Sections |. 0. and [II]. 






V. Enss 


Inverse Scattering 


7 Time Scales for Interaction and Spreading 


For high energy states as constructed in (|3.8| ) scattering theory becomes simple because 
two time scales, a short interaction time T/(jp) and a longer kinematical time of spread¬ 
ing T Sp { p) satisfy T/(p)/T Sp (p) -> 0 as |p| —>► oo. For a potential which is essentially 
supported in a ball of radius R the interaction time is of the order T/(p) = R/v( p). 


Due to propagation like ( |3.11| ) we have 


S(t + , t_) w S if both ±t ± »T/(p). (7.1) 

More precisely, for Tp, $1 as in ( |3.8| ) and any e > 0 there is a radius p(e) such that 
uniformly for large |p| (which satisfy ( |3.10|) ) 


|(4' [5 — S(t + ,t_) ] 4p)| < 


u(p) 


if ± f± > p(s)/v(p). 


(7.2) 


Intuitively, p(e) measures the radius of the interaction region and the extension in 


configuration space of the states up to effects of size e. The bound (|7.2|) is physically 
intuitive and it is a crucial estimate which will be used in Section ^ to justify the 
interchanging of limits. 

To prove ( fh2|) one has to bound a term 


lim ('■ ^c it Ho [r/o+u s ] ^ h^Hq 

p p 


' dt ± e it[H 0 +V] e -itH 0 & 

dt p 


< / dt V s e 


s —it Ho (f./ 


«pI 


(7.3) 


and a similar term with t_ and Tp . The integrand in( |7.3| ) can be split into two terms 
\\V S e ~ itH ° $'-11 


< IV' F(|x| > !>(p)f/2)H + ||V*|| ||F(|x| > v(p)t/2) e-“«° 42 


=: fti(o(p) t) + h 2 { v(p) t) 

where F(-) here denotes the multiplication operator with the characteristic function 
of the indicated region in configuration space. The functions h\ and h 2 are integrable 
due to (|4.2|) and (|3.11|) , respectively. With the new variable r := u(p) t the integral 
(|F3|) is bounded by 


v{P) 


dr [h\(r) + h 2 {r ) ] < 


u(p) 


for r + > p{e) large enough. This proves ( |F2|) . 


The kinematical time scale of spreading Ts p ( p) denotes the time after which 
spreading of wave packets becomes relevant in the time evolution. As 

tfo(p) Tp = #o(p) e ?:px T 0 = e ipx H 0 ( p + p) T 0 

we will expand the kinetic energy function around p 

H 0 (p + P) =: H 0 (p) + Vi7 0 (p) ■ P + H 2 (p, p). (7.4) 
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The first summand is a number giving an irrelevant phase, the second equals v(p) • p 
by (|3.9| ). It is the dominant term which - as a multiple of the momentum operator 
- generates a translation of the wave packet without changing its shape. Only the 
third term H 2 (which is defined by ( [f.4|) ) is responsible for the spreading of the wave 
packet. In our examples of “power like” Hamiltonians this part of the Hamiltonian is 
weak compared to the translational component: On a compact subset of momentum 
space like p E supp © 0 

T r(P) \H 2 (P, p)| ^ const n frj ^ 

T Sp { p) u(p) |p| IpI-oo 

In the nonrelativistic case we have H- 2 (p } p) = p 2 /2 m which is independent of p and 
bounded on supp 0o • 

|#2(p, p)| = p 2 /2m = p 2 /2 
u(p) |p| /m |p| 

In the relativistic case the denominator is bounded but the numerator decreases. With 
the shorthand a := \/p 2 + m 2 we get for |p| large enough 

\H 2 (P, p)| = |a{\/l + 2 (pp/a 2 ) + (p 2 /a 2 ) - 1 - pp/a 2 }| < 0 p 2 /2 
v(p) |p| /a ~ |p| 

This proves dZD- Therefore, the time T$p( p) is by a factor proportional to |p| longer 
than T/( p). For large |p| we may choose times when the scattering due to the potential 
is over but the spreading has not yet really started. Alternatively, we may use radii for 
this splitting like the interaction radius p(e) of (|7.2|). 

Usually, an interacting time evolution is complicated because the translation of a 
wave packet, its spreading, and the influence of the potential all occur at the same time. 
In the high-energy limit it is sufficient for the calculation of the scattering operator to 
treat translation of wave packets rather than their correct free evolution. Since in this 
limit spreading occurs only when the interaction is negligible, i.e. when the free and 
interacting time evolutions are almost the same, the effect of spreading is cancelled 
(becomes invisible) in the scattering operator. Thus, high energy scattering is simple 
and it can be inverted simply! The motion during the interaction is dominated by 
the translational part which is common to classical and quantum physics. The typical 
quantum effect of spreading of wave functions which results from the absence of localized 
states with sharp momentum is of lower order. 

8 High Energy Scattering 

The crucial uniformity of the estimate ( |7. 2|) enables us to interchange the limits 
±t± —> 00 and |p| —> 00 . This simplifies the remaining discussion very much. Actually, 
as to be expected, not the time but the separation from the region of a strong potential 
determines the quality of approximation. With correspondingly chosen variables 
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r± t±v( p) we have 


lim ($1, S$p) = lim lim (& p , S{t + ,t_)$ p ) 


IpHoo 


= lim lim $1, S 

|p|->oo =tr±— >o ' D 


= lim lim $L 5 

±r±—>oo |p|—>oo 


|p|—>oo ±t±—>oo' P 


n(p) ’ u(p) 


v(p) ’ v(p) 






( 8 . 1 ) 


As seen in ( |7.2| ) the asymptotic equality ( |8.1| ) remains true even after multiplication 
with u(p) which is a much stronger statement in the nonrelativistic case. To determine 


(Op, p) = (*i, e-^S^t 

for large finite times and p consider e.g. the second pair of factors in 

e -«px e it_[H 0 +V s ] e ~it_H 0 e *px 
_ e it_[i?o(p+p)+^ s (x)] e ~it_H 0 (p+p) 

— [Ho(p)+v(p)-p+H 2 (p,p)+V s (x)] e ~it_ [H 0 (p)+v(p)-p+H 2 (p,p)] 

— [v(p)-p+H 2 (p,p)+V s (x)] e ~it_ [v(p)-p+H 2 (p,p)] 

_ e ir_ [a;-p+{i/ 2 (p,p)/i)(p)}+{V s (x)/i;(p)}] g -ir_ [oj p+{ff 2 (p,p)/u(p)}] 


( 8 . 2 ) 


(8.3) 


using again t± = r±/v( p) and the direction w = v(p)/u(p) as in ( |3.9|) . Due to (|7.5| ) 
the functions of the momentum operator 


[w-p + {if 2 (p,p)Mp)}] —-P (8.4) 

|p|—>oo 

converge in strong resolvent sense and similarly for the other exponent. Therefore, for 
fixed r_ and large |p| the following approximation is good: 

e ir_ [u:-p+{H2 (p,p)/i)(p)}+{U s (x)/t!(p)}] g -ir_ [u>-p+{.ff 2 (p,p)/i;(p)}] 


e %r_ [oj-p+{V s (x)/i)(p)}] -ir_ui -p 


(8.5) 


= exp 


—i 


v(p) 


dr H s (x + us r) 


( 8 . 6 ) 


The approximation (|8.5|) is the only approximation we have to make! If {i/ 2 (p, p)/u(p)} 
would commute with (H s (x)/T(p)} then we would have exact cancellation and (|8.5|) 
would be an equality as well. A careful estimate of the correction terms can be given 
for all Hamiltonians considered here. It is uniform in r and when compared to 
"(V s (x)/u(p)} it has additional falloff like l/|p| for p — 2 00 due to (0)- For the 
proofs we refer to the papers cited in the theorems below. 

Equation (|3.6|) is verified easily because as functions of r_ both expressions solve 
the same initial value problem 
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The same analysis of the first two factors in the expression (|5.2|) for the scattering 
operator yields analogously to (|8.6[) the factor 


exp 


dr E s (x + co r ) 


TIPJ Jo 

Combining this with ( |8.6| ) we obtain for large |p| 

e^S' r r 


;VP X 


$o 


$o > eX P 


V(p) 


w(p)' v(p), 

dr E s (x + ur) } $ 0 


(8.7) 


(%, S*p)« K, exp 


V{P) J -c 


dr V s (x + cor) > $ 0 ) . 


( 8 . 8 ) 


These expressions confirm our intuitive expectation discussed in Section |3] that the lim¬ 
iting behavior of the scattering operator depends strongly on the growth or boundedness 
of u(p) as |p| —> oo. We will discuss the two cases separately. 


9 High Energy Limits of the Scattering Operator, 
the Short-Range Case 

Next we give the limiting behavior of the scattering operator in several typical cases, 
$p, & p , and per as given in (|3.8|) , (|3.10|) . The integrals extend over the real line. In 
the quotations we include similar results obtained by other methods, sometimes under 
more restrictive assumptions. 


Theorem 9.1 (scalar relativistic, short-range, ||14| ) 

For the scalar relativistic Hamiltonian 

H = y/p 2 + m 2 + V s (x) 
with u(p) —> 1 one obtains 

(t'p, S J>p) (%. exp j-i J dr V‘(x + wr) j i’« ) . 


(9.1) 


If, however, u(p) —> oo we can expand the exponential in (|8.8|) 


exp 


v(p) 


dr V s (x + cor) 


1 - 


u(pj 


dr V s (x + cor) + 


(9.2) 


which explains the following nonrelativistic result. The leading behavior of the scatter¬ 
ing operator is the identity operator (no scattering). The next order correction depends 
on the potential. 
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Theorem 9.2 (nonrelativistic, short-range, [|10 

For the Hamiltonian 

ff =^ p2 + v "M 


18 


v(p) ($p, i{S - 1) $p) 


($' 0 , r(x + wr)$ 0 ). 


(9.3) 


This result is to be expected from the Born approximation. It holds also under the 
given weaker assumptions on the falloff of the potential where the validity of the Born 
approximation is not established. 

The estimate (|7.2|) and the remark following ( |8.6|) justify that multiplication with 
u(p) ~ |p| is permitted. The approximation is better than C s /u(p) . 


Remark 

In all these limits there are error bounds for large but finite |p| which are explicit. 
E.g. in equation ( |9.3| ) (and similarly in (|10.5|) etc.) we obtain 


«(P) (%, KS- l)*p) 


I dr ($' 0 , E s (x + cor) <h 0 ) 


const ($ 0 , $ 0 , V s ) 
|p| 


Theorem 9.3 (Dirac equation, [|12|1 , ||14j| ) 

For the Dirac Hamiltonian of a particle in a continuous short-range electromagnetic 
field 


H — ap + j3m + V(x), V(x) = (A 0 — ol- A)(x) 
the high-energy limit of the scattering operator is 





dr (A 0 — ct uj 


cj-A)(x + cor) 



In the Fouldy-Wouthuysen representation let S± denote the scattering operator on the 
positive/negative energy subspaces and x NW the Newton- Wigner position operator then 


s-lim e * S±e l p * N w 

|p|—>oo 


exp 


-i / dr (A 0 la'W oj- A)(x w + wr) 


Similar results hold for other matrix valued potentials V. 


10 High Energy Limits of the Scattering Operator, 
Long-Range Potentials, Magnetic Fields, etc. 


The only long-range potential in physics is the Coulomb potential between electri¬ 
cally charged particles. The mathematical treatment includes more general potentials 
which may decrease slower and need not be centrally symmetric. A convenient class 
of long-range potentials are twice continuously differentiable functions with C £ (x) —> 
0 as Ixl —► cx) and 


1 / 



VC £ (x)| < const(1 + |x|)-( 3 / 2 )- 
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|AE £ (x)| < const (1 + |x|)- 2 " £ , e > 0j. (10.1) 

Local singularities may be treated with the short-range part. If long-range potentials 


are present the ordinary scattering operator (|5.2|) , ( |5.3|) no longer exits, the free time 
evolution has to be replaced by a better approximation generated by H 0 (t ) := H 0 + 
V\tp/m). The factor e~ lt - H ° in ( fn^) has to be replaced by 


U D (t_, 0) := exp — it_ H 0 — i I dtV e (tp/m) 


( 10 . 2 ) 


and similarly U D (0, t + ) for the term with t + . With this Dollard correction the modihed 
scattering operator exists: 


S D (t + ,t_) := [U D (0,t + )T e^ H ° +v '+ vl ' e «-[W+^] f jz>( t _,0) , 


(10.3) 


S := s-lim S (t ,t_), satisfies S' 

t . —>oo 


(10.4) 


t —>■—oo 


The splitting of the potential into its short- and long-range parts is not unique. Different 
choices - as long as (fL2|) and (|10.1|) are satisfied - correspond to different labellings of 
asymptotic configurations and are physically equivalent. 

To define a modihed scattering operator some knowledge of the long-range behavior 
of the potential is needed. A typical physical situation is that the total charge is known 
but the charge distribution is to be determined by scattering experiments. 

Theorem 10.1 (nonrelativistic, long-range, []8|) 

For a splitting of the potential as given above let the Hamiltonian be 

H = i p2 + l/,(x)+ 

and let S D be the corresponding Dollard modified scattering operator then 


u(p) ($1, i(S D - 1) <Lp) - dr ($(,, [ W(x + wr) - V\un) ] <h 0 ) 


♦ / dr ($(,, W(x + wr) $ 0 ). 


IpHoo 


(10.5) 


Note that the integral f dr [ V f (x + ujr) — V\ujr)] converges for every x. Another 
version of the theorem yields at once the full potential by using the high-energy limit 
of [pi, S D ] for any modihed scattering operator. 

If one has a magnetic held described by a vector potential A in a nonrelativistic 
Schrodinger operator then the leading term in the interaction is A- p/m ~ A-v(p). 
Inserting this into (|8.8|) suggests that one obtains a nontrivial limit of the scattering 
operator although the kinematics is nonrelativistic. The other potentials show up as a 
correction to the leading behavior of lower order. 
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Theorem 10.2 (nonrelativistic with magnetic fields and potentials of short 
and long range, |l5|, 0, [|]) 

Let the vector potential A(x) be a short-range multiplication operator and 
H = ^-[ p-A(x)] 2 + F s (x) + ^(x) 

( 10 . 6 ) 

= -^p 2 — A — + ^-{i div A + |A| 2 ) + 

2 m m 2m 

The scattering operator itself has the limit 

dr co ■ A(x + cor) 


s-lim e * px Se ipx — eX p 

|p|—>oo 


(10.7) 


which yields a unique magnetic field F(x) = rfA(x) for continuous A(x), and - after 
fixing a gauge - the vector potential A(x). In a second step one obtains as a lower 
order correction to the leading term 




jGj g —i f dr a>-A(xH-u>r) _ 


$p ) — known terms(A, V e , $o, $ 0 ) 


|p|—>oo 


> / dr (<I> 0 , F s (x + u;r) $ 0 ). 


( 10 . 8 ) 


Clearly, only the magnetic field gLA(x) (= curl A(x) inz/ = 3 dimensions) is a physical 
quantity. However, fixing a gauge and choosing a vector potential is needed to write 
down the Schrodinger operator (|10.6|) (choice of a representation). For homogeneous 
magnetic fields see a, a. 

Let us now turn to nonrelativistic multiparticle scattering systems with short- and 
long-range potentials. As the interaction becomes weak for high speeds the incoming 
and outgoing scattering channels are the same in the high-energy limit. Excitations 
and rearrangements are of lower order. 


Theorem 10.3 (nonrelativistic multiparticle systems) 

In two-cluster scattering the limit of the channel scattering operator for high relative 
speed of the clusters yields the effective intercluster potential between the two bounded 
subsystems, ||. 

In the totally free channel select a pair of particles with high relative speed and all 
other particles far away (or with even higher speeds). Then the pair potential of this 
pair is obtained from the limit. Multiparticle potentials can be obtained iteratively after 
the pair potentials are known, @ , 0 . 


The method has been extended and applied to other problems as well, see e.g. 

time-dependent IV-body Schrodinger operators [^1 
and nonlinear Schrodinger equations 


Stark effect Hamiltonians [20 


time-dependent Dirac operators 


11 Reconstruction of the Potential 

The condition v > 2 {multi dimensional inverse problem) enters here to obtain from 
the above limits reconstruction formulas and uniqueness. For bounded continuous (or 
more general) functions V s the expression 

W s (x,co) \= f drh s (x + wr) (11.1) 
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is the X-ray transform of V s . In v — 2 dimensions lines and hyperplanes are the same. 
Therefore, (|11.1|) is the Radon transform as well. The latter is known to be uniquely 
invertible because the assumption (|4.2| ) implies V s € L 2 (M 2 ), see e.g. Theorem 2.17 
in Chapter I of [|ll|]. The inverse Radon transform yields the unique potential. In 
higher dimensions one fixes e.g. X 3 , ..., x v and reconstructs the “slices” subsequently. 
In particular, it is sufficient to vary lo in a two dimensional plane. For unbounded 


or discontinuous potentials the expectation value between states from a dense set of 
nice vectors (like those which satisfy ( |3.7|) ) effectively smoothes the potential. This is 
enough to reconstruct the potential as a multiplication operator, see e.g. ||. 

The reconstruction of the magnetic field is similar but more complicated, see |IJ, 
1111 for details. 


or 
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